Inhomogeneous quantum diffusion and decay of a meta-stable state 
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We consider the quantum stochastic dynamics of a system whose interaction with the reservoir 
is considered to be linear in bath co-ordinates but nonlinear in system co-ordinates. The role of 
the space-dependent friction and diffusion has been examined in the decay rate of a particle from a 
meta-stable well. We show how the decay rate can be hindered by inhomogeneous dissipation due 
nonlinear system-bath coupling strength. 



I. INTRODUCTION 

The system-reservoir modelpl H H 0, IE 0] has 

been the standard paradigm for description of quantum 
dissipation over the last few decades. Its application is 
ubiquitous in physics and chemistry. In the overwhelm- 
ing majority of the situations the interaction between 
the system and reservoir is considered to be linear in 
bath co-ordinates and linear in system co-ordinates. This 
in turn relates the additive noise of the thermal bath 
with linear dissipation of the system through fluctuation- 
dissipation relation. However when coupling between the 
reservoir is linear in bath co-ordinates but nonlinear in 
system coordinates , on e encounters multiplicative noise 
U H 113 IH 13 111 El El 13 113 and nonlinear dissipa- 
tion in the form of co-ordinate dependent friction. The 
role of the space-dependent friction in classical context 
has been explored in several issues, e. g., charge transfer 
reaction in a polar medium and activated rate processes 
in the overdamped regime 

[3 13 13 El E3 fluctuation- 
induced transport |23l l24l . stochastic resonance prj |27|. 
noise-induced transitioned l2flj . etc. to name a few. 

In the present paper we are concerned with space- 
dependent friction and multiplicative noise in a quantum 
mechanical context. A standard treatment of quantum 
dissipation based on linear interaction between the sys- 
tem and the reservoir wasput forward in early eighties by 
Caldeira and Leggett H 0] , which found extensions and 
wide applications in several areas of condensed matter 
and chemical physics HHHI3] particularly in quantum 
interference devices, polaron in a magnetic field, quan- 
tum decay processes etc. However, the exploration of 
physics of nonlinear system-reservoir coupling in a quan- 
tum system has been the subject of recent interest only. 
It has been observed that quantum dissipation can re- 
duce the appearance of meta-stable state and barrier 
drift in a double- well potential[2l|. Nonlinear coupling 
has been modeled by Tanimura and co-workers |5j in 
the study of elastic and inelastic relaxation mechanisms, 
and their interplay in the Raman and infrared spectra. 
A quantum transport as a consequence state-dependent 
diffusion has been investigated|l7j in a periodic poten- 
tial. A closer look into the nature of nonlinear coupling 
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(say, a linear in bath co-ordinate but quadratic in system 
co-ordinate) in contrast to linear coupling suggests that 
while the latter induces displacement (Fig. 1(a)), the for- 
mer coupling induces fluctuations of the system poten- 
tial (Fig. 1(b)). The effect of this fluctuation has been 
addressed in both s pin- spin and spin-lattice relaxation 
processes |30l l3lt 133. l33 | . The fluctuation in the system 
potential due to nonlinear coupling, in general, tends to 
make the effective potential nonlocal in addition to mak- 
ing the diffusive process inhomogeneous in nature. This 
consideration is therefore of significant relevance in the 
context of quantum rate processes and it is worthwhile to 
examine the nature of the decay of a meta-stable state as 
a result of inhomogeneous quantum diffusion. Based on 
a coherent state representation of the noise operator and 
canonical thermal Wigner distribution 34] of the bath os- 
cillators nonlinearly coupled to the system we have re- 
cently developed the quantum stochastic dynamics in the 
overdamped limit |17|. This has been carried out after 
adiabatic elimination of the momentum in the relevant 
Langevin dynamics which takes care of systematic quan- 
tum corrections to the potential. 

Our aim in this paper is to search for the quantum 
signature of this nonlinear coupling on the decay rate, 
particularly in the pre-exponential factor and the gener- 
alized effective potential. We have emphasized how to 
take care of correct order of dissipation and appropriate 
thermodynamic consistency to avoid spurious term in the 
steady state distribution^, HH] and show how the decay 
rate can be reduced by inhomogeneous dissipation. 

The outlay of the paper is as follows: In Sec. II we 
introduce briefly the basic aspects of quantum stochastic 
dynamics with nonlinear coupling on a general setting. 
This is followed by a quantum Langevin equation in the 
overdamped limit. Sec. Ill is devoted to the calculation of 
decay rate from a meta-stable state under the influence 
of inhomogeneous diffusion and dissipation. This paper 
is concluded in Sec. IV. 



II. QUANTUM STOCHASTIC DYNAMICS 
WITH NON-LINEAR COUPLING 



We consider a particle of unit mass coupled to a 
medium comprised of a set of harmonic oscillators with 
frequency ujj . This is described by the following system- 
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bath Hamiltonian 



H = 



(2.1) 



Here q and p are the coordinate and momentum op- 
erators of the particle and the {xj,pj} are the set of 
coordinate and momentum operators for the bath oscil- 
lators with unit mass. The system particle is coupled 
to the bath oscillators nonlincarly through the general 
coupling terms ;j- /(<?)■ Cj is the coupling strength. The 

classical counterpart [35| of the form H2.lt is known for 
many years and also the nonlinear coupling in quantum 
system has been studied in a few occasions 0, . The 
potential V(q) is due to external force field for the sys- 
tem particle. The coordinate and momentum operators 
follow the usual commutation relations [q,p\ = ih and 
[xj,pk] = ihSjk- 

Eliminating the bath degrees of freedom in the usual 
way we obtain the operator Langevin equation 0, Il7| for 
the particle 



q(t) = p(t) 



(2.2) 



p(t) = -V'(q(t)) - /'(?(*)) / /W)b(* ~ t')p(t')dt 



(2.3) 



-fwmt) 



where the noise operator fj{t) and the memory kernel 7(f) 
are given by 



m = 



J2 \ — %(0) - f(q(0)) cosu 3 t 



-Pj(0) sin uijt 



and 



lit) 



(2.4) 



(2.5) 



It is clear from the operator Langevin equation 
Ea. l|2.3(l for the system that the noise operator is mul- 
tiplicative and the dissipative term is nonlinear with re- 
spect to system coordinate due to the nonlinear coupling 
term in the system-bath Hamiltonian. In the case of lin- 
ear coupling, i.e., f(q) — q Ea. H2.3fl reduces to a quantum 
generalized Langevin equation |7] in which the noise term 
is additive and the dissipative term is linear. 

In the Markovian limit the generalized quantum 
Langevin equation Ea. (|2.3ll reduces to the form 



?(*) - Pit) (2.6a) 
P{t) = -V(q(t))-rlf(q(t))} 2 p(t) 

+rmmt) (2.6b) 



where L is dissipation constant in the Markovian limit. 

Following Ref.@ we then carry out a quantum mechan- 
ical average (..) over the product separable bath modes 
with coherent states and system mode with an arbitrary 
state at t = in Ea. H2.6al) and Ea. (|2.6bl) to obtain a 
generalized Langevin equation as 



q = p 



p = -V'(q)+Q v -r[f'(q)] 2 p + Q 1 
+f'(q)v(t) + Q2 



(2.7a) 



(2.7b) 



where the quantum mechanical mean value of posi- 
tion and momentum operators (q) = q and (p) = p. 
Qv — V'{q) — (V'(q)) represent quantum correction 
due to nonlinearity of the system potential. Q\ — 
T [[f'iq)?P ([f'(q)?P)]_ and Q 2 = rj(t) [</'($)) f'(q)} 
are also quantum correction terms due to nonlinearity of 
the system-bath coupling function. 

Furthermore quantum mechanical mean Langevin 
force is given by, 



v(t) 



2 12 

^ (%(0))-(/(g(0)))l^cos^ 

C 3 ) U j 



(2. 



To realize rj(t) as an effective c-number noise we now 
introduce the ansatz |6j, l34l 1371 1381 139L |40j that the mo- 

mentum (pj(0)) and the shifted coordinates {— (xj(0)) — 

{f{q(0)))} of the bath oscillators are distributed accord- 
ing to a Wigner canonical thermal distribution of Gaus- 
sian form as: 



^[{^(x.m-ifimm^m)^ 

l <K(0)) 2 + 4{4 <* 3 (0)M/(«(0))>} 2 ] 



Afexp 



2huj (nj {uij )+ ^ ) 



(2.9) 



so that for any quantum mechanical mean value, 

Oj ({f (%(0)> - (/(9(0)))}, (^(0))) which is a func- 

tion of mean value of the bath operators (xj(0)) and 
(pj(Q)}, its statistical average (...)s is 



(Oi) g = / OjVj d(f>M) d{^{xM) - </($(o))» 

J c j 

Here rij(u>j) indicates the average thermal photon 
number of the j-th oscillator at the temperature T and 
fij(ujj) — [exp(-^y) — arL d A/" is the normalization 
constant. 

The distribution Vj (Ea. H2.9fl ) and the definition of 
statistical average Eg . (f2. lUf) imply that c-number noise 
r)(t) must satisfy 



(V(t))s = 



(2.11) 



3 



(V(t)v(t'))s = % ( coth 2^r) C0S ^ ( * " ^ 



In the Markovian limit the noise correlation becomesl2z 

lai 



(V(t)ri(f))s = 2D S(t-t') 
D 1 



Thw I n(uo) + ^ 



(2.13a) 
(2.13b) 



where luq is the average bath frequency and the spectral 
density function is considered in the Ohmic limit. The 
use of Gaussian form [Ea. l|2.9|l ] for a bath distribution 
needs some elaboration. It follows from a simple calcu- 
lation of partition function of a nonlinearly coupled har- 
monic oscillators, that a Gaussian form of ansatz is not 
valid, in general. In the present case we emphasize that 
one is concerned here with a system-reservoir nonlinear 
coupling where the coupling is nonlinear in system coor- 
dinates but still linear in bath coordinates, so that the 
use of coherent state basis for harmonic oscillators enable 
us to factor out the dependence of nonlinear system co- 
ordinate in our calculation. This is also reflected in the 
fluctuation-dissipation relation (|2.12() which retains its 
usual form, in spite of nonlinear system-reservoir interac- 
tion. Further to this we emphasize that the overdamped 
situation in our treatment refers to large L rather than 
r [f'{x)} 2 . For a better validity of the theory the nonlin- 
earity should not be too strong. 

The equations (|2.11l) and (|2.12l) imply that the c- 
number noise rj(t) is such that it is zero centered and sat- 
isfies the standard fluctuation-dissipation relation. We 
thus identify rj(t) as a classical looking noise with quan- 
tum mechanical content. The quantum Langevin equa- 
tion can now be rewritten as follows: 



q = p 



(2.14) 



V = -V'(q) + Q v -T[f(q)] 2 p-2Tpf(q)Q f -TpQ 3 
- 2Tf(q)Q 4 - T Q 5 + f'(q) V (t) + Q f V (t) (2.15) 

The above equation is characterized by a classical force 
term, V , as well as its correction Qv ■ The terms con- 
taining r are nonlinear dissipative terms where Qf, Q3, 
Q4 and Q§ are due to associated quantum contribution in 
addition to classical nonlinear dissipative term r[/' '(q)] 2 p. 
The explicit expressions for Qv, Qf, Q3, Qa and Q5 are 
given in the Appendix-A. The last term in the above 
equation contains a quantum multiplicative noise term in 
addition to the usual classical contribution f'{q)rj(t). The 
classical limit of the above equation was derived earlier 
by Lindenberg and Seshadri |35j. Furthermore quantum 
dispersions due to potential and coupling terms in the 
Hamiltonian are entangled with nonlinearity. The quan- 
tum noise due heat bath on the other hand is expressed 
in terms of the fluctuation-dissipation relation. 



While considering hydrodynamic interaction, i.e., 
when the fluctuation is position/state dependent or 
equivalently when the noise is multiplicative with respect 
to system variables the conventional adiabatic reduction 
of fast variables does not work correctly To obtain a 
correct equilibrium distribution Sancho et al had pro- 
posed an alternative approach to Langevin equation in 
the case of multiplicative noise system. By carrying out 
a systematic expansion of the relevant variables in pow- 
ers of r _1 and neglecting terms smaller than 0(r _1 ) they 
obtained a Langevin equation corresponding to a Fokker- 
Planck equation in position space. This description leads 
to the correct stationary probability distribution of the 
system with coordinate-dependent friction. 

We follow the same procedure in the quantum mechan- 
ical context. In this limit the transient correction terms 
Qi and Q5 do not affect the dynamics of the position 
which varies in a much slower time scale. So the quan- 
tum Langevin equation Ea. (|2.14|l and Eq. l|2.15|) can be 
written as, respectively, 



q = p 

p = -V'(q)+> 



) V -T h(q)p + g{q) V (t) 



where 



Hq) = [f'(q)} 2 + 2f(q)Q f 
9(q) = f(q)+Qf 



(2.16) 
(2.17) 



(2.18) 
(2.19) 



The variable g(q) arises due to nonlinearity of the system- 
bath coupling function f(q), Qf being the quantum cor- 
rection to the classical contribution /'((?). For a linear 
coupling function g(q) reduces to a constant. 

Following the method of Sancho et al j| further we 
obtain the Fokker-Planck equation in position space cor- 
responding to Langevin equation Ea. (|2.17|l 



dP(q,t) d 



Of 



V'(q) ^ Qv 



+ ^°7T 
oq 

+ D l 



Th(q) 
1 



T (h(q)) 2 g ' {q)9{q \ 
g(q) d g{q) 



r h(q) dq T h{q) 



P(q,t) 
P(q,t) 



(2.20) 



The Stratonovich prescription leads to corresponding 
Langevin equation as given by 



V'(q) ~ Qv 
Th(q) 



D 



g'(q)g(q) , g(q) 



r r h( q ) 



v(t) (2.21) 



Eq. I|2.21|) is quantum Langevin equation for multiplica- 
tive noise with position dependent friction in the over- 
damped limit (i.e., corrected upto 0(l/r)). 

In the classical limit, i.e., haj <C ksT, h(q) = [f'{q)\ 2 , 
g(q) = f'(q), Qv = and Dq = T ksT so the quantum 
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Langevin equation reduces to its classical counterpart as 
derived by Sancho et al 9] 



1 



q = 



in. 



r [f'(qW 



-V'(q) - r k B T 



f"(q) 
f'(q) 



/'(?) #) 

(2.22) 



DECAY RATE FROM A METASTABLE 
STATE 



Attention is now restricted here to the problem of 
barrier crossing dynamics in the presence of nonlinear 
system-bath coupling. The particle coordinate q (which 
in our case is the quantum mechanical mean position) 
corresponds to the reaction coordinate and its values 
at the minima of V(q) denote the reactant and prod- 
uct states separated by a finite barrier, the top being a 
metastable state representing the transition state. We 
choose the potential and the coupling function of the fol- 
lowing forms ; 



V(q) = l -aq 2 - l -bq z 
/(<?) = h q + i& 2 q 2 



(3.1) 
(3.2) 



a and b are the potential parameters and the maximum 
and the minimum of the potential are at the point q = 
qb = a/b and q = q a = 0, respectively. To clarify the 
role of the nonlinear system-bath coupling we now turn 
to Ea. (|2.1(l . in which system-bath interaction is given by 
(ujj is absorbed in Cj), 



Hint — — 



^2^3 x 3 f{q) 



= -{c\Xi + c 2 x 2 + ...)(b 1 q+ -b 2 q 2 ) (3.3) 

The linear-linear coupling term Cj b\q) is proportional 
to b\ while square- linear coupling term (^ Cj b 2 q 2 ) is 
proportional to b 2 . The linear-linear coupling model a 
very commonly used mechanism for studying quantum 
dissipative dynamics, e. g., energy dissipation from vi- 
brational system mode to the heat bath modes during 
population decay. In addition to energy dissipation due 
to linear-linear coupling, one also encounters fluctuation 
in the system potential due to nonlinear coupling. We 
refer to Fig. 1(b) where a schematic illustration of the 
fluctuation of system-bath coupling in a relevant poten- 
tial well is presented. Since this fluctuation is the key 
mechanism for state-dependent diffusion, it is apparent 
that the calculation of transition rate in presence of state- 
dependent diffusion differs from the calculation of escape 
rate for particles subject to thermally uniform noise, in 
that we may conveniently substitute the role of potential 
V(q) by a generalized potential ^(g) corresponding to 
the former case. A simple way to determine transition 
rate is to consider a steady state condition. The Fokker- 
Planck equation Ea. H2.20[) in the overdamped limit can 



be rewritten in a more compact form as 



dP{q,t) d 1 



dt 



dq r %) 



V'(q) ~Q V + 



Do d g{qf 



r 8q h{q) 



P(q,t) 
(3.4) 



In the overdamped limit the stationary current from 
Ea. (|3.4|) is given by 



J = — 



r%) 



VI \ n d (^) 2 \ 

V[q) - Qv + -Tq{W)) 



Pst(q) 
(3.5) 



The Eq. H3.5f) can be rearranged into the following form, 



— <{ P st L(q) exp 



*(g) 
D n 



"TT %) ex P 



*(g) 



where 
*(<?) 



(3.6) 



q V'W) ~ Qv 



dq'; 



and D q = ^ . To calculate a steady state solution we as- 
sume an absorbing boundary at a coordinate q s past the 
intervening potential maximum at q — So integrating 
Eq. (|3.t)|) between q = q a and q = q s , we have 



P s t{q a ) 


L(q a ) exp 


r*(9a)i 

D q J 


r<ls r 


h(q) exp 




dq 



(3.8) 



The steady escape rate over the top of the potential is 
given by flux over population, 



n„ 



(3.9) 



where n a is the population around the potential mini- 
mum (q a )- The zero current defines an equilibrium con- 
dition in the neighborhood of q a . Thus with J = 0, prob- 
ability distribution function near the potential minimum 
is given by 



Pst(q)=P s t(q a )^y4exv 

L (q) 



(3.10) 

In terms of the generalized potential Ea. H3.10() can be 
rewritten as 



-*(ff) , *(g«) 



D a 



(3.11) 



Population around the bottom of the potential is given 
by 



Pst(Q) dq 



Pst{0a) L{q a ) exp 



\*(q a y 


M2 l 


[-*(?)] 


[ D q _ 







(3.12) 
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where qi and g 2 are two points around potential mini- 
mum. So the expression for transition rate is given by 



1 



1 



r In" d 1 Hi) ex P 



~*(g) 



(3.13) 

By linearizing the generalized potential around maximum 
and minimum and evaluating the integrals in the usual 
way we obtain 



|tt"( ga )|* \*"(q b )\?L(q a ) 
k = M%) 6XP 

can be simplified further to yield 
q V'(q')-Q v 



D„ 



(3.14) 



*(«/) 



dq' (V'(q') - Qv) 1 + 



(W) + Q/) a ; 

(3.15) 



The rate constant l|3.14l) then reduces to the following 
expression, 



k = 



V[^(l + Ai)-A„ 1 ][^(l + A 2 )-A„ a ] 



x exp 



27rr h(q b ) (1 + Ai 

A, 



(3.16) 



where w a = V"(g a ) and w fc = V"{lb) are the frequencies 
at the bottom and barrier top, respectively. A 1; A 2 , A V1 
and A V2 are the contributions due to quantum correction 
and are given by, 



Ai = 



Q 3 - Q) 
(/'(?) + Qff 



A 2 = 



Q 3 - gf 

(/'(<z) + Q/) : 



J 9=9!, 

(3.17) 



& V1 = 



A,. 



-1 




-J 




dq) 


Qy ( 



(/'(«) + Q f f 

Q 3 -Q} 
(f'(q) + Qf) 2 



g=Qa 



(3.18) 



The quantum nature of the transition rate Ij3.16|l is man- 
ifested through the quantum correction due to nonlinear- 
ity of the potential and coupling function and the quan- 
tum inhomogeneous diffusion. It also contains the signa- 
ture of nonlinear coupling in the pre-exponential factor 
and in the generalized potential ^(q). More generally, 
the effect of state-dependent diffusion makes its presence 
felt in the dynamics of decay. 

In order to check some inherent consistency of the rate 
expression we now examine the following situations. 



(a) First consider that the coupling function is linear, 
i. e., b 2 = and f(q) = q (bi is assumed to be unity). 
Then the generalized potential is given by 



9(q) 



dq(V'(q') - Qv) 



(3.19) 



and the transition rate reduces to 



x exp 



27iT 
V(q b ) - V{q a ) 



exp 



Acj) 
~D~„ 



where 



Acj) 



in, 



Qv(q) 
D n 



dq 



(3.20) 



(3.21) 



This quantify the characteristic quantum decay rate from 
a meta-stable state for particles subject to thermally uni- 
form noise. On the other hand as the classical limit is 
approached all the quantum correction terms (A Vl , A V2 
and A(f>) due to nonlincarity of the potential tend to van- 
ish and the quantum coefficient (D q ) gets replaced by the 
classical one {D q — > k b T). The transition rate (|3.2U|) re- 
duces exactly to Kramers rate in the overdamped limit. 



k = ^T eXP 



V(q a ) 



k h T 



(3.22) 



(b) In order to check the thermodynamic consistency 
we consider further a symmetric periodic potential with 
periodicity 2ir, i.e., V(q) = V(q+2ir) and periodic deriva- 
tive of coupling function with the same periodicity as that 
of the potential, i.e., f'{q) — /'(g + 27r). Such a potential 
had been subject of interest in the problem of classical [23j 
and quantum transport |l7j| in inhomogeneous media. 

Since the potential is periodic, Qv is also a periodic 
function as Qv = V(q) — (V(q)). Similarly Qf is also 
periodic as Q f — (f'(q)) -/'(g), and also Q 3 + 2f'(q)Qj 
is also periodic since Q 3 +2f(q)Q f = {[f(q)] 2 )-[f'(q)} 2 . 
So from Eas. (|2.18|l and (|2.19[) it is clear that h{q) and g(q) 
are also periodic functions of q with periodicity 27r. Now 
it is easy to check that *ff(q) is also a periodic function 
with periodicity 2tt. 



$(q) = W( 9 + 2tt) 



(3.23) 



Consider the bottom of the potential well at q a = 2im 
and potential maxima at q b = (2n + l)n. Then with the 
help of Eq. I|3.16[) we find a transition rate from valley n 
to valley n + 1 



vn-\-l,n 



V[^(l + Ai)-A U1 K(1 + A 2 )-A 



x exp 



2irT h(q b ) (1 + Ai) 
D n 



(3.24) 
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Repeating these considerations for the transition rate 
from valley n + 1 to valley n yields 



+ 1 



VK(1 + A 1 )-A i;i ][^(l + A 2 )-A„ a ] 



x exp 



2irTh{q b ) (1 + Ai) 

gXgb) - ^(gq +2?r) 



(3.25) 



Thus the average particle velocity to a preferential direc- 
tion is given by 



,dq 

{-rr} = 2?T ( fcjj+i „ — k n n+ i J 



27rfc n+ i,„ ( 1 - exp 




*(g a ) - (g Q + 2vr) 



(3.26) 



We thus conclude that there is no occurrence of current 
for a periodic potential and periodic derivative of cou- 
pling with same periodicity. At the macroscopic level this 
confirms that there is no generation of perpetual motion 
of second kind, i.e., no violation of second law of ther- 
modynamics. Therefore the thermodynamic consistency 
based on symmetry considerations ensures the validity of 
the present formalism. 

Before we proceed to numerical results it is important 
to clarify the issues regarding the calculation of quan- 
tum correction terms. The details of the calculations of 
quantum correction terms are shown in the Appendix-A. 
One can calculate the value of quantum dispersion terms 
(Sq n ) by direct numerical simulation of the coupled Eqs. 
(A10) subject to appropriate boundary conditions. It 
is also instructive to deal with the quantum correction 
terms in the analytical way to find out the approximate 
value of quantum dispersion terms. For overdamped limit 
we neglect the 5p term from Eq. i|A9|l to obtain 



Sq 



1 



dt T [f'{q)f 

1 



r [f'W 



[-V"(q) 5q-2Tpf(q)f"( q )S q } 
W)f"{q)5q] + 0(5?) (3.27) 



With the help of Eq. (|3.27|l we then obtain the equa- 
tions for (5q n ) 



r lf'(q)} 2 
4 

2 



[-V"{q){5f)} 

[Tpf'{q)f"{q){5q 2 )} 

[ V (t)f'(q)(6q 2 )]+0({Sq 3 }) 



(3.28) 



dt [q} r[/'(g)]3 
6 



[-V"(q){Sf)] 

[I>/W(?)(<5<f}] 



r [f'(qW 

r - 9 Ht)f"(q)(Sq 3 )]+0({6q 4 )) 



r [f'(qW 



(3.29) 



and so on. (It is apparent from Eqs. i|3.28!3.29|) that in 
the overdamped limit the higher order quantum contri- 
butions are small since each successive order of correction 
is lower than the preceding one by a factor of -jt) 

A simplified expression for the leading order quan- 
tum correction term (5q 2 ) can be estimated by neglecting 
the higher order coupling terms in the square bracket in 
Eq. (|3.28|l (since the nonlinearity of the potential is small 
the terms of the order of f"(x) are small) and rewrit- 
ing it as d(Sq 2 ) = r Tjr^m V"(q)(Sq 2 } dt. The over- 
damped deterministic mean motion on the other hand 
gives dq = — r rprop ^ These together yield after 
integration^ 



(5q 2 )=A q [V'( q )f 



(3.30) 



where A n 



" [^'(go)V anc ^ ^° * S a Q^ 11 ^ 111 mechanical 
mean position at which (Sq 2 ) becomes minimum, i. e., 
(Sq 2 ) qa — Tjh/ujQ, ujo being defined earlier. 

The relevant quantum correction terms IjAlj! , l|A4jl and 
(|A5|) can be rewritten for the potential i|3.1[l and coupling 
term (|3.2|l as follows, 



1 



Q f 
Q 3 



V"'(q)(Sq 2 ) 







= [/"(<z)]W> 



(3.31) 
(3.32) 
(3.33) 



V""(q) and f"'(q) vanish for the chosen metastable po- 
tential and square-linear coupling. So by calculating the 
quantum dispersion term {Sq 2 ) we can find out the values 
of Ai, A 2 , A V1 , A V2 and A</> as given by Eas (|3.17l3.18|) 
and <pOT|) . The corrections (|3~3T|) and lpH3")l are reminis- 
cent of the perturbation corrrections to Kramers escape 
rate due to color It may also be mentioned that im- 
portant nonlinear terms are those represented by a linear 
chain coupled to a linear heat bath[44| and in spite of the 
different coupling Hamiltonian, the effects on the relevant 
rates are similar. 

We now proceed to illustrate numerically the behavior 
of quantum transition rate given by Eq. (|3.16|l . It appears 
that in addition to the barrier height the pre-factor is 
also affected by the quantum correction terms. The cal- 
culated quantum transition rate Ea. l3.16|) implies that 
both the barrier height and the frequency factor contain 
the effects of quantum corrections due to nonlinearity 
of system potential and coupling function. The effect of 
quantization of the reservoir and quantum correction due 
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to system nonlinearity is apparent in Fig. 2 in the vari- 
ation of In A; with l/T (Arrhenius plot) for varied con- 
tribution of nonlinear coupling. All the curves exhibit 
linearity at higher temperature and nonlinear variation 
is observed at lower temperature. Decrease in hi (square- 
linear coupling contribution is proportional to 62) results 
in enhancement of transition rate. This physically im- 
plies that with increased values of hi potential well be- 
comes comparatively steeper (as shown in the Fig.l) and 
hence the transition rate decreases. 

The effect of quantization of a classical transition rate 
is shown in Fig. 3, where we make a comparison of k 
vs l/T profiles for classical and quantum cases. One 
observes that at low temperature regime the classical 
transition rate is significantly lower in magnitude than 
quantum rate and at higher temperature the effect of 
quantization become insignificant. To clarify the effects 
of quantization we present also an Arrhenius plot com- 
paring the classical and quantum regimes in the inset of 
Fig. 3. From the inset plot it is clear that at the high tem- 
perature regime the plot exhibits linearity in the both 
cases, which is the standard Arrhenius classical result. 
In low temperature limit, however, one observes a much 
slower variation in the quantum case than that of clas- 
sical one. To examine this in more detail it is not dif- 
ficult to see that the exponential factor in the rate ex- 
pression Ea. H3.16|) can be reduced to a form containing 

the usual Arrhenius term exp — f tlSkL^tiSz) ^ times a 
T 2 enhancement term at low temperature of the form 



exp 



fax> 



( ^(q b )-Mq a ) \ 
\ 2(fc7T ) 



The full quantum behavior can 

be interpreted in terms of an interplay between the quan- 
tum diffusion coefficient D q and the quantum correction 
due to nonlinearity of the system potential appearing in 
^(g). When the temperature of the system is very low, 
i. e., in the vacuum limit or in the deep tunneling region 
the anharmonic terms in the potential do not contribute 
significantly. On the other hand as the temperature of 
the system increases significantly, D q increases resulting 
in decrease of the effective potential and hence D q and Q 
(quantum correction terms) compete to cancel the effect 
of each other at higher temperature. 



IV. CONCLUSION 

Based a theory of stochastic dynamics of a quantum 
particle in inhomogeneous media under overdamped con- 
dition we have calculated the rate of decay of a meta- 
stable state of the system under the influence of quan- 
tum state dependent friction and diffusion. It has been 
shown that the nonlinear interaction between the sys- 
tem and bath have its imprints on the pre-exponential 
factor as well as on the form of a generalized effective 
potential. Furthermore the nonlinearity of the effective 
potential can be identified as a typical consequence of 
quantum effect. The present theory can also be extended 
to other selected areas of chemical physics where the 



fluctuations in the potential in the system-reservoir cou- 
pling are important. From a purely quantum mechanical 
point of view such interactions give rise to multi-photon 
transition terms in the dynamics. Thus their contribu- 
tion to phase relaxation in addition to population decay 
is of considerable significance. The difference between 
linear-linear and quadratic-linear coupling in the system- 
reservoir interactions terms is likely to be also important 
if one calculates the time evolution of flux-flux correla- 
tion functions |45|| which in consequence gives rise to time- 
dependent rate constant. The focus of the present arti- 
cle being the non-equilibrium steady state rate constant 
and ensuring a correct description of equilibrium in the 
quantum system which often remains vague in many re- 
lated situations, we feel that an elaboration on this issue 
needs more detailed consideration. We note, in passing, 
that the treatment presented here concerns stochastic 
processes under overdamped condition in the Markovian 
limit. Extension of the theory to non-Markovian and 
weak friction regime is worth-pursuing in the context of 
rate theory and related contexts. 
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APPENDIX A: QUANTUM CORRECTION 
TERMS 



Referring to the quantum nature of the system in the 
Heisenberg picture we now write the system operators q 
and p as q = q + Sq and p — p + Sp respectively. Sq and Sp 
represent quantum fluctuations around their respective 
mean values. By construction {Sq) = (Sp) = and they 
also follow the usual commutation relation [Sq, Sp] = ih. 
Using a Taylor series expansion in Sq around q and Sp 
around p we express Qv, Qi, Q2 and Q3, Qi, Q5, and 
Qf as functions of q, p, (Sq n ) and (Sp n ) as follows; 

Qv = V'(q) - (V'(q)) 

= -E^ n+1 w<^"> ( A1 ) 

n>2 

Q 1 = T[[f\q)fp-([f\q)Yp)] 

= -T[2pf(q)Q f + P Q 3 + 2f'(q)Q i + Q^A2) 
Q2 = v(t) [(/'(?)> -/'(«)] 

= V(t) Qf (A3) 
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where, 
Qf 

Qs 



n>2 

^^i T l/™ +1 (< z )/"+ 1 ( (Z )(^ m ^") (A5) 
L — ' L — ' ml n\ 

m. > 1 n > 1 

Eir +1 (?)W"^) (A6) 



an infinite hierarchy of equations. Up to third order 
we construct, for example, the following set of equations 
(A a) which are coupled to quantum Langevin equations from 
1 ' EH ; 



E E —Af m+1 (l)f n+1 (q)(Sq m Sq n Sp) 



(A7) 



The dynamics of these correction terms can be 
calculated 01 with the help of the following equations, 
which can be derived using the operator Langevin equa- 
tions H2.6afl and 12.6b|l and by carrying out quantum me- 
chanical average over the initial product separable coher- 
ent bath states. 

5q = Sp (A8) 

dp = -V"{q)5q - ]T -.V n+1 (q) [Sq n - (Sq n )} 
c — ' rv. 

- 2Tf( q )f"( q )Sq + 2Tf'{q) £ -J n+1 (q) (Sq n ) 



n>2 



dt 



(5q 2 ) = (5q5p + 5p6q) 



(AlOa) 



dt 



(SqSp + SpSq) = -2 X (q,p)(Sq 2 ) + 2(Sq 2 



> 2 > 



T{f'(q)} 2 (5q5p + 6p8q)-aq,p)(6q 3 ) 
2Tf(q)f'(q)(6q 2 5p + SpSq 2 ) (AlOb) 

-2T[f( q )} 2 (dp 2 ) - X (q, P) (SqSp + Sp6q) 



dt 



dt 



- \({q,p)(Sq 2 Sp + SpSq 2 ) 

- 2Tf'( q )f"(q)(Sq6p 2 + 6p 2 5q) 

(Sf) = I(6q 2 6p + Sp5q 2 ) 

(6f) = -3r[f( q )} 2 (dp 3 ) 



:X(q,P){SqSp 2 + Sp 2 5 q ) 



(AlOc) 
(AlOd) 

(AlOe) 



dt 



(6q 2 5p + SpSq 2 ) 



-2 X (q,p)(5f) + 2(6q5p 2 + 6p 2 5q) 
r[f( q )] 2 (6q 2 Sp + 6p6q 2 ) (AlOf) 



E E —<-J m+ \q)P + \q) [^q m Sq n - (Sq m Sq n )]p T M5f + Sp 2 Sq) = 2(Sf) - 4 X (q,p)(Sq 2 Sp + SpSq 2 ) 

- 2T[f'(q)} 2 (SqSp 2 +Sp 2 Sq) (AlOg) 



v[f'(q)] 2 5p + 2f{ q ) Y, -J n+ \q) [Sq n sp - (Sq n sp)} 



dt 
where 



n>l 



E E — } -J m+1 (q)f n+1 (q) {6 q m 6q n Sp-(6 q m 6q n Sp)} 



i>l n>l 



to! n 



f"(q)Sq+y2^f n+1 (qM n -{Sq n ) 



(A9) 



X(q,p) = V"( q )+2Tpf( q )f"( q )-r 1 (t)f"( q ) 

(AlOh) 

C(9,P) = V"'( q ) + 2Tpf'( q )f>"( q ) 

+ 2Tp[f'( q )} 2 - v (t)f"'( q ) (AlOi) 



The operator correction equations can be used to yield For other details we refer to • 
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Figure Captions 

Fig.l: Schematic illustration of the effect of linear- 
linear(a) and square-lincar(b) system-bath coupling on 
a relevant potential system. 

Fig. 2: Plot of In A; vs. 1/T for different strength of 
nonlinear coupling & 2 ((i)&2 = 0.01 (dashed line) (ii) b 2 = 
0.1 (solid line) (iii) b 2 = 0.2 (dotted line) for a fixed 
parameter set a = 0.5, b — 0.15, b\ = 1.0 and T = 1.0 
(all the quantities are in dimensionless unit). 

Fig. 3: A comparison between classical(solid line) and 
quantum (dotted line) decay rate plotting k as a func- 
tion of 1/T for the parameter set a = 0.5, b = 0.15, b\ = 
1.0, b 2 = 0.10 and T = 1.0. The inset presents a varia- 
tion of In k with 1 /T comparing classical (solid line) and 
quantum (dotted line) cases for the same parameter set 
as in the main figure, (all the quantities are in dimension- 
less unit). 
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